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Abstract 

A liquid Argon TPC (ICARUS-like) has the ability to detect clean neutrino 
bursts fi'om type-II supernova collapses. In this paper, we consider for the first 
time the four possible detectable channels, namely, the elastic scattering on elec- 
trons from all neutrino species, Vf, charged current absorption on Ar with pro- 
duction of excited v^, charged current absorption on Ar with production of 
excited CI and neutral current interactions on Ar from all neutrino flavors. We 
compute the total rates and energy spectra of supernova neutrino events includ- 
ing the effects of the three-flavor neutrino oscillation with matter effects in the 
propagation in the supernova. Results show a dramatic dependence on the os- 
cillation parameters and in the energy spectrum, especially for charged-current 
events. The shock breakout phase has also been investigated using recent simu- 
lations of the core collapse supernova. We stress the importance of the neutral 
current signal to decouple supernova from neutrino oscillation physics. 
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1 Introduction 



Core collapse supernovae are a huge source of all flavor neutrinos. Neutrino astrophysics 
entered a new phase with the detection of neutrinos from the supernova SN1987A in 
the Large Magellanic Cloud by the Kamiokande and 1MB detectors [1]. In spite of 
these fundamental neutrino observations, the 19 events observed are not statistically 
significant enough to obtain precise quantitative information on the neutrino spectrum. 
Currently running neutrino detectors like Superkamiokande or SNO have the capabil- 
ities to provide high statistics information about supernova and neutrino properties if 
a supernova collapse were to take place in the near future. 

The flavor composition, energy spectrum and time structure of the neutrino burst 
from a galactic supernova can give information about the explosion mechanism and 
the mechanisms of proto neutron star cooling. In addition, because neutrinos arrive 
before any other signal from a supernova, it is possible to provide an early alert to the 
astronomical community combining the observations of several experiments. This will 
allow an early observation of the flrst stages of the supernova explosion. 

The neutrino signal from a galactic supernova can also give information about the 
intrinsic properties of the neutrino such as flavor oscillations. Although new data from 
solar, atmospheric, reactor and accelerator neutrinos [2] have contributed to the under- 
standing of the neutrino properties, still the neutrino mixing angle dis and the nature 
of the mass hierarchy, i.e., normal or inverted, remain unknown. These parameters can 
be probed by the observation of supernovae neutrino bursts, since neutrinos will travel 
long distances before reaching the Earth and will, as they travel through the exploding 
star, in addition traverse regions of different matter densities where matter enhanced 
oscillations will take place. In the case of the 613 angle, matter enhancement has the 
striking feature that very small mixing angles, beyond any value detectable by next 
generation accelerators, could in fact alter signiflcantly the neutrino spectrum. Hence, 
supernova neutrinos provide indeed a complementary tool to study the 613 angle. 

The main question of course is to understand to which extend can the supernova 
and the neutrino physics be decoupled in the observation of a single supernova? On 
one hand, the understanding of the supernova explosion mechanism is still plagued by 
uncertainties which have an impact on the precision with which one can predict time, 
energy and flavor-dependent neutrino fluxes. On the other hand, the intrinsic neutrino 
properties are not fully known, since the type of mass hierarchy and the value of the 
6'i3 angle are unknown, and in fact large uncertainty still exist on the prediction of the 
actual effect of neutrino oscillation. 

There have been various studies on future supernova neutrino detection and their 
interplay with neutrino oscillations physics. Dighe and Smirnov [3] estimated qualita- 
tively the effects of neutrino oscillations in a collapse-driven supernova on the neutron- 

ization peak, the distortion of energy spectra and the Earth matter effect. Dutta et 
al [4] and also Takahashi et al [5] showed numerically that the total number of events 
for Superkamiokande and SNO detectors increases dramatically when there is neutrino 
mixing. 
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Detailed supernova simulations (see e.g. [6]) indicate that the average neutrino en- 
ergy is flavor dependent and hierarchical (E,^^ < Ep^ < E^-^^^p^^). Great uncertainties 
still exist in the spectral and temporal evolution of the neutrino fluxes. Although the 
hierarchy of the average energies of the z/g, z/g and p^, is believed to hold, quan- 
titatively their specific spectra remain a matter of detailed calculations (see [7], and 
[8] and references therein). Under the assumption of hierarchical neutrinos, neutrino 
flavor oscillations should make the spectra of and Pe harder (see [9] and references 
therein). This will lead to a change in the neutrino observations in the detector because 
the cross sections depend strongly on the neutrino energy and they are different for 
different flavors. 

It is hence well known that uncertainties in the supernova explosion mechanism 
and the effect of neutrino oscillations can change dramatically the number of observed 
events in a given experiment. Indeed, this dependence can be used to determine the 
supernova and mixing parameters, once the next supernova has been detected. Most 
detailed studies in the literature however concern water or heavy-water targets (see 
e.g. [10, 11, 12, 13]). In addition, some authors advocate the use of speciflc variables 
in order to disentangle supernova from neutrino oscillation physics (see e.g. [9]). 

In this work we study the capabilities of a liquid Argon TPC (ICARUS-like) to 
detect neutrinos from supernova collapses and its possible contributions to this held. 
Our study is motivated by the proposed 3 kton ICARUS detector [14] that should 
become operational at the Gran Sasso Underground Laboratory and by ideas for a 
very large liquid Argon TPC with a mass in the range of 100 kton's (see e.g. [15]). 

We think that a liquid Argon detector has excellent potentialities for supernova 
physics, namely because of its high granularity and bubble-chamber like features, all 
channels, namely elastic scattering off atomic electrons, as well as charged and neutral 
nuclear scattering processes are detectable and identiflable. The simultaneous observa- 
tion of nuclear charged and neutral currents provide high-statistics and golden channels 
to decouple supernova from neutrino oscillation physics. 

It is known that precise elastic scattering can also yield information on the flavor 
content of the flux (see e.g. [2] in the context of the solar neutrino puzzle), because of 
different cross-sections depending on ffavor and helicity. However, this requires good 
statistical precision that is more difficult to achieve with target electrons than with 
nuclear target. Hence, nuclear processes are expected to contribute most. 

A Liquid Argon TPC can also separate the various channels by a classification 

of the associated photons from the K, CI or Ar de-excitation, which exhibit specific 
spectral lines, or by the absence of associated photons in case of the elastic scattering. 
It is hence possible to detect z/g, z/g and other flavors independently. 

In section 2, we briefly recall the general framework we used in this study. In par- 
ticular, we identify flve scenario for the cooling of the supernova. Neutrino oscillations 
in vacuum and in the supernova matter are discussed in section 3. 

In section 4, we discuss the relevant neutrino cross-section on Argon and identify 
the experimentally detectable classes of events. New calculations on the cross section 
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for the absorption processes based in the random phase approximation (RPA) were 
used allowing computing for the first time the contribution of the Pg '^^Ar absorption 
reaction and of neutral current reactions. 

In section 5, a particular interest is taken on detection of the early phase neutrino 
signal from supernova in liquid Argon. We investigate the expected events coming from 
the first hundreds of milliseconds after the supernova bounce using the calculations on 
the shock breakout process in core collapse supernovae and taking neutrino oscillations 
into account. 

In section 6, we compute the expected neutrino rates in liquid Argon, for different 
primary supernova neutrino fluxes scenarios. Moreover, the energy spectra modified by 
oscillations are plotted for different conditions on the oscillation parameters. From su- 
pernova neutrino observations it is possible to extract information about the character 
of the mass hierarchy and the ^13 mixing angle. 

2 General framework and supernova explosion sce- 
narios 

When the iron core of a massive star (M > 8 M©) overcomes its hydrodynamical 
stability limit (the Chandrasckhar mass), the core collapses raising its density up to 
many times the nuclear density. During this infall stage of the collapse, a first i^e 
burst is emitted ("m/a// hursf), since the high density of matter enhances the electron 
capture by protons. The total energy radiated during this phase is roughly 10^^ ergs. 

This emission does not continue indefinitely, since at a density ~ 10^^ g cm~^ 
neutrinos are trapped in the stellar core and go to equilibrium with matter. The 
anomalous density produces an elastic bounce of the core, which results in a shock 
wave. After the core bounce^ also neutrinos of other flavors begin to be produced. The 
neutrinos are trapped in a region ( "neutrinosphere" ) whose size is different for different 
neutrino flavors. A deeper neutrinosphere corresponds to a higher temperature. 

The wave propagates through the star and loses energy in dissociating nucleons. 
When the shock crosses the Vf, neutrinosphere, an intense burst of v^. ( "shock breakouf 
or "neutronization bursf) is produced by electron capture on the large number of 
protons liberated by the shock. The characteristic time of the breakout burst is 3- 
10 ms and the total energy radiated in i/g neutrinos during breakout is ~ 3 x 10^^ 
ergs. 

After the breakout stage, the luminosity rapidly decreases, while the luminosities 
of other flavor neutrinos increase. At the end of the neutrino diffusion inside the mantle 
the energy is practically equafly distributed between the various neutrino flavors. The 
shock wave stalls without driving off the stellar mantle and envelope, z/g and Ug are 
absorbed on the nucleons liberated by the shock. Such processes supply new energy to 
the wave, which is revived ~ 500 ms after the bounce. This energy transfer is known 
as "neutrino heating" . Continued mass accretion and convection below the shock wave 
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also deposit energy ans thus contribute to the shock revival. The reinforced shock can 
propagate within the stellar matter and expel the external layers into the space. 

After the explosion the star loses energy, mainly by neutrino emission, and cools 
down, forming a neutron star or a black hole. This final stage {^^coolinff^) requires ~ 
10 s and takes away > 99 % of the gravitational energy of the star (~ 2-4 x 10^^ ergs). 

In order to obtain quantitative estimates, the following general considerations are 
taken into account for the present analysis: 

• SN type-II at a distance of 10 kpc (galactic supernovae) 

• Gravitational binding energy of 3 x 10^^ ergs (radiated away as neutrinos). 

In our analysis of the event rates and energy spectra, we simplify the supernova neutrino 
emission process and assume that it can be split into two distinct phases: the early Ue 
pulse expected in the first 40 ms of the collapse {shock breakout) and the rest beyond 
the first 40 ms which we treat as a cooling stage, which could in reality last up to ~ 
10 s. 

• For the shock breakout we use the time-dependent energy spectra from Burrows 

et al. The calculation is based on dynamical models of core collapse supernovae in 
one spatial dimension, employing a newly-developed Boltzmann neutrino radia- 
tion transport algorithm, coupled to Newtonian Lagrangean hydrodynamics and 
a consistent high-density nuclear equation of state. Full details of the simulation 
can be found in [16]. Recent studies on the initial mass progenitor dependence 
seem to indicate that the major features of the early neutrino burst should be 
almost independent of the initial mass of the star [8] . 

• For the cooling stage wc assume that the energy spectra of neutrinos is a Fermi- 
Dirac distribution. In order to take into accoimt the uncertainties in the mod- 
elling of the cooling phase, we consider five different scenarios which correspond 
to simulations found in the literature. The various scenarios can be parameter- 
ized by different assumptions on the average energies and relative luminosities of 
the different neutrino flavors. Table 1 summarizes the different possibilities for 
the primary supernova neutrino flux parameters and neutrino average energies 
considered in this work. 

We assume the following neutrino oscillation situations: 

• Normal (mf < m| << m|) and inverted (m| << mf < m|) mass hierarchies are 
investigated in the analysis. 

• Matter effects inside the supernova are included. 

• The following neutrino oscillations parameters are used in the analysis, obtained 
from other neutrino observations [2]: 

sin^ 2^23 = 1- Amis ^ Amii = 3 x lO'^eV^ 

sin^ 012 = 0.3 Amii = 7 x lO-^e^^ (1) 

sin^ ^13 < 0.02 
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Cooling 
scenario 


MeV 


{E.J 
MeV 


MeV 


Luminosity 


Reference 


I 


11 


16 


25 




[17] 


II 


13 


16 


23 


Lj/e = ^Pe = ^Ux 


[10] 


III 


13 


16 


17.6 




[7, 18] 


IV 


13 


16 


17.6 


Li/e = Lpe = 2Li/^ 


[7, 18] 


V 


13 


16 


17.6 


Li/g = Lp^ = 0.5 Lj,^ 


[7, 18] 



Table 1: Supernova cooling phase scenarios. 



3 Neutrino oscillations and matter effects 

Neutrinos produced in the high density region of the iron core interact with matter 
before emerging from the supernova. Due to the non-zero masses and non-zero vac- 
uum mixing angles among various neutrino flavors, flavor conversions can occur in 
supernovae. 

The transitions are produced mainly in the so-called resonance layers. The reso- 
nance density can be computed as [3]: 



TTlN 

coszO (2) 



2V2Gf E Ye 



where Gp is the Fermi constant, Am^ is the mass squared difference, 9 is the mixing 
angle, E is the neutrino energy, uin is the nucleon mass and Y^. is the number of 
electrons per nucleon. 

Since the inner supernova core is too dense to allow resonance conversion, we can 
consider two resonance points in the outer supernova envelope: one at high density 
p Ki 10^-10^ g cm~^ (H-resonance) which is governed by the atmospheric parameters 
(Am|]^ and ^13) and the other one at low density p 10-30 g cm~^ (L-resonance) , 
characterized by the solar parameters (Am|i and ^12). 

The transitions in the two resonance layers can be considered independently and 

each transition is reduced to a two neutrino oscillation problem. The H-resonance 
lies in the neutrino channel for normal mass hierarchy and in the antineutrino channel 
for the inverted hierarchy. The L-resonance lies in the neutrino channel for both the 
hierarchies [3]. 

The probability of transition between one neutrino eigenstate to another in the 
resonance layer (L or H) is called jump probability (Pl or Ph)- 
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Figure 1: Jump probability in the H-resonance, Ph, as a function of sin^ ^13 for three 
values of the neutrino energy (left), and as a function of the neutrino energy for different 
values of sin^ 6'i3 (right). 

The propagation through the low density region is always adiabatic for the LMA 
solar parameters. It means that Pl and Pl are both equal to and the neutrino mass 
eigenstates remain unchanged. 

Considering a matter density profile of the star p(r) oc r^^, the jump probability in 
the H-resonance can be written as [3]: 



Ph oc exp 



-const sin^ 9i ' '^^ 



as 



(3) 



Indeed, the neutrino conversion in the high density region depends on the mixing angle 
6'i3 and the mass squared difference between the involved flavors, as shown in figure 
1. The variation of the jump probability with the 13-mixing angle and the neutrino 
energy is plotted. We see that the crossing probability increases with energy. Three 
regions can be distinguished: 

1. Non adiabatic region (sin^6'i3 < 2 x 10~®): For these values of 613 the jump 
probability is almost equal to 1, independently of the neutrino energy. 

2. Intermediate region (2 X 10-6 < sin2 0i3 < 3 X 10^^): In this region the value 
of Ph, between and 1, depends on 6*13 and neutrino energy as expressed in the 
equation 3. 

3. Adiabatic region (sin^ 6^13 > 3 x 10^^): The jump probability is equal to for any 
value of the neutrino energy. Then, there is no conversion between eigenstates in 
the resonance layer. 
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The effect of oscillations in the neutrino fluxes is summarized in equation 4, where 
the expected fluxes at the detector {(f)^) are computed in terms of the original fluxes 
at the SN core (0°) and the survival probabilities Pee — P{^e ^e) and Pee = P(i^e 



0.. = C:Pee + C(l-^ 



ee ) 
ee ) 



+ d,^ = dfjl - Pee) + 0° (1 + Pee) 



eel 



with 0°^ = (f)°^ — (f)l^ — = = If the original fluxes were equal between 
flavors, any oscillation effect could be observed. 

The survival probabilities can be written in terms of jump probabilities as [3] : 

(5) 

L\'-^e'2\ 

for normal hierarchy (Am32 > 0) and 



Pee = PHPLjUell' + P^ (1 " ^l) | ^e2 P + (1 " P^jl^esP 
Pee=(l-PL)|^7eip + Pi|[/e2p 



Pee^PL\Uel\' + il-PLm 



e1\ 



|2 



Pee = P^/(l - PLWelf + PLPi/|C/e2P + (1 - PhW. 



e3 



(6) 



for inverted hierarchy (Am|2 < 0). 

Considering the values of the jump probabilities and using the standard parametriza- 
tion of the mixing matrix, the survival probabilities can be computed as shown in figure 
2. The evolution of the probability Pee (Pee) that a neutrino emitted as v^. (v^ at the 
neutrinosphere remains as v^. (z/e) at Earth is plotted as a function of the ^13 mixing 
angle for different values of neutrino energy. 

The upper plot corresponds to the normal hierarchy case (n.h.) and the lower one 
to the inverted hierarchy (i.h.). For n.h., the survival probability in the antineutrino 
channel does not depend on the neutrino energy and only very weakly on Qyi- How- 
ever, the neutrino channel is very sensitive to ^13 via the Pu parameter. For inverted 
hierarchy the neutrino channel is independent on ^13 while the antineutrino channel 
depends on the 13-mixing angle and neutrino energy. 

According to the adiabaticity conditions in the H resonance and the value of the 
6*13 angle, three regions can be distinguished. If energy spectra are different, then 
and the mass hierarchy can be probed as follows: 

• For small mixing angle (sin^ ^13 < 2 x 10~^), there are no effects on ^13 and both 
hierarchies give similar results. Therefore, in this region we can not distinguish 
among hierarchies and only an upper bound on sin^ ^13 can be set. 

• For intermediate ^13 (2 x 10~^ < sin^ ^13 < 3 x 10~^) maximal sensitivity to the 
angle is achieved and measurements of the angle are possible in this region. 



^We neglect the possible effect of the propagation through the mass of the Earth. 
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Figure 2: Survival probabilities for normal (top) and inverted (bottom) hierarchies 
a function of sin^ ^13 and different values of neutrino energies. 
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• For large mixing angle (sin^ 6'i3 > 3 x 10"^) maximal conversions occur. The effect 
of the mixing is strong in neutrino (antineutrino) channel if the mass hierarchy 
is normal (inverted). Wc arc able to probe the mass hierarchy in this region but 
only a lower bound on 6*13 can be put. 



4 Cross-sections and expected supernova neutrino 
signal on liquid Ar detectors 

At low energy, where nuclear processes dominate, three nuclei-scattering processes can 
be studied (z/g charged current absorption on Ar with production of excited K, v^. 
charged current absorption on Ar with production of excited CI and neutral current 
interactions on Ar from all neutrino flavors) in addition to elastic scattering on atomic 
electrons from all neutrino species. There is a possibility to separate the various chan- 
nels by a classification of the associated photons from the K, CI or Ar de-excitation, 
which exhibit specific spectral lines, or by the absence of associated photons in case of 
the elastic scattering. 

In a hquid Argon TPC, it is hence possible to identify the following events: 



1. elastic scattering on electrons 



(-) _ (-) _ /„x 

V e ^ V e (J) 



2. chctrged- current (CC) interactions on eirgon 

i/e ^°>lr ^ e- "^K* (8) 
z7, ^Mr ^ e+ (9) 

3. neutral-current (NC) interactions on cirgon 

V' 40^^ 40^^* ^IQ^ 



4.1 Elastic scattering processes 

As well known, the elastic neutrino scattering off target-electrons has a total cross 
section that increases linearly with energy: 

a{vee- Uee-) = 9.20 x 10-*^E^^{MeV) cia 

a{vee-^Vee ) = 3.83 x IQ-^^Ep^jMeV) cm^ , . 

(T{i^t.,re- ^iy^,,re-) = 1.57 x IQ-^^^.^^^^MeV) cm' ^ ^ 

<^{i^^^,re- ^ iy^,re-) = 1.29xlO-^^E^^^{MeV) cm' 
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All neutrino species contribute to elastic scattering. The experimental signature 
consists in a single recoil electron. Since, the direction of this electron is highly cor- 
related to the incoming v direction, these events have the potentiality of precisely 
determining the location of the supernova source[19]. 

Figure 3 shows the cross sections of the elastic processes as a function of neutrino 
energy. We recall that the reaction involving z/^'s is the most probable since both CC 
and NC currents contribute. Electron antineutrinos are suppressed by helicity. The 
cross-section for muon and tau neutrinos is the smallest, since they proceed through 
NC current only. 

These features are important since while the nature of the neutrino cannot be 
identified in elastic events, the total number of expected elastic events for a given flux 
depends on their nature because of its cross-section dependence. Hence, if the total 
flux of neutrino is constrained, the number of elastic events is sensitive to the nature 
of the neutrinos emitted by the supernova. 

4.2 Nuclear processes 

In order to compute the event rates, the relevant nuclear cross-sections must be in 
principle known with relatively good precision from energies ranging from a few MeV 
up to 100 MeV. For the nuclear processes, we use the calculations computed for the 
first time using the RPA including all the multipoles [20]. 

Figure 3 shows the cross sections of all the processes as a function of neutrino 
energy. The NC processes present high cross sections of the same order of magnitude 
as the ones from the CC reactions. They are even bigger than those from the PgCC 
events. This can be understood in terms of the energy threshold and nuclear effects of 
the various reactions. The Q-value of the CC is 1.50 MeV while for the CC case 
is 7.48 MeV, hence electron antineutrinos are kinematically disfavored. In addition, 
the l"*" and 0"*" allowed transitions are suppressed by nuclear structure reasons and the 
cross section is dominated by the other multipoles. Hence, electron antineutrinos cross- 
section is smaller than for neutrinos. The first excited state of Argon has an energy 
of 1.46 MeV with respect to the ground state, so the threshold for neutral current 
reaction is low. 

In the case of the CC process, various estimates are available in the literature 
(see e.g. [21] and [22]). It was assumed in the past, that the vJ^^Ar absorption rate 
proceeds through two main channels: the superallowed Fermi transition to the 4.38 
MeV excited isobaric analog ^°K* state, and other Gamow- Teller transitions to several 
excited states. This assumption is correct for neutrino energies relevant to solar 
neutrino detection, however is not adequate for the energy range of supernova neutrinos, 
which extends up to 100 MeV. In the recent calculations that we use, the cross-section 
has been computed by random phase approximation (RPA) for neutrino energies up to 
100 MeV, including all the multipoles up to J=6 and both parities. Figure 4 compares 
this cross section from three different sources: the dashed line corresponds to the shell 
model calculations performed by Ormand [21] considering only Fermi and Gamow- 
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Neutrino energy (MeV) 

Figure 3: Neutrino cross sections relevant to the supernovae detection with a hquid 
Argon TPC detector. 
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Teller transitions; the dotted line is the assumption made in [16] where they take the 
total absorption cross section as 3 times the cross section of the pure Fermi transition 
computed in [22]; finally the solid fine is the cross section used in this analysis [20]. 

We see that the assumption of atot = 3 itf overestimates the cross section for energies 
smaller than 30 MeV, which is the energy range of z/g from breakout, and therefore 
affects the expected number of events from the burst (discussed in section 5). 

On the other hand, both Ormand and atot = 3 cross sections underestimate 
the cross section for high energies up to 100 MeV, relevant for supernova neutrinos 
from the cooling phase (discussed in section 6). The cross section including all the 
multipoles grows faster than the others, showing that indeed not only the Fermi and 
Gamow-Teller transitions are important but also the rest of multipoles. 

5 Supernova collapse and neutrino breakout phase 

5.1 Time and energy dependence without neutrino oscilla- 
tions 

The neutrino shock breakout burst is the signal event in the supernova core collapse 
evolution. Its detection and characterization could test fundamental aspects of the cur- 
rent collapse supernova paradigm. The prompt electron neutrino burst is in principle 
observable and represents a diagnostic of the fundamental collapse supernova behavior. 
The characteristic time of the breakout burst is of the order of tens of ms. 

The advantage of the early phase analysis is that it should not be affected by the 
time evolution of the density structure of the star nor whether the remnant is a neutron 
star or a black hole (see e.g. [16]). The shock wave takes about 2 seconds to reach the H- 
resonance region. Therefore, the potential time dependence of neutrino oscillations due 
to shock propagation can be neglected in this phase. Clean detection of the early part 
of the supernova signal is therefore a very important aspect of experimental supernova 
detection. 

As already mentioned, a liquid Argon TPC has an excellent sensitivity to neu- 
trinos mainly through the CC process. This provide unique information about the 
early breakout pulse. We will however show that neutral currents as well as elastic 
scattering, although lesser in rate, can provide important information as well on the 
shock breakout. 

In order to study the detection of the shock breakout, we compute the differential 
number of neutrinos detected at a distance D from the supernova: 

m.,.,. = ^(^yE.),E.M (12) 

where Ntargets is the number of targets, D is the supernova distance, Ei, is the neutrino 
energy, t is the detector time, dN,^/ dE^dt is the time and energy spectra of the neutrinos 
and (t{Ei,) is the cross section of a reaction process. 
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Figure 4: CC cross section as a function of the neutrino energy. The dashed hne 
corresponds to the Ormand [21] cross section calculation, dotted line assumes that the 
total cross section of the absorption interaction is 3 times the cross section of the Fermi 
transition [16] and the solid line is the cross section used in this analysis calculated 
from RPA including all the transitions [20]. 
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We use the simulation of the early phase of the neutrino burst performed by Burrows 
et al. [16] which yields (fN^,/ dE^^dt in discrete steps of time. 

The prediction for the time (energy integrated) evolution of the neutrino fluxes 
for the different flavors in absence of neutrino oscillations is shown in Figure 5 (top). 
The time origin is defined such that the supernova core bounce time corresponds to 
t = 200 ms. This time is taken as the reference for the rest of the processes. The 
Figure 5 (bottom-left) shows the time-integrated energy spectrum of the first 470 ms of 
the time, namely up to 270 ms after the bounce. As expected, the flux is characterized 
by the bright peak of the emission due to the supernova shock breakout. 

For comparison, we show in Figure 5(bottom-right) the particular energy distribu- 
tion at 100 ms after the bounce {t = 300 ms). This picture illustrates how the energy 
spectra of the different neutrino flavors is predicted to change with time and how the 
bright z/g becomes superseded by the neutrinos of all other flavors and helicities. 

For the sake of calculations, we deflne the shock breakout burst by integrating up 
to t — 240 ms, namely we define a duration for the burst of 40 ms. We use these in 
the next section. In order to simplify the analysis, we treat the times after t — 240 ms 
as part of the coohng phase (see Section 6) . 

5.2 Effect of flavor oscillations and predicted event rates 

Neutrino oscillations modify dramatically the fiavor content of the neutrinos emitted 
from the supernova during the burst. A reduction of the peak due to the effect of os- 
cillations can be readily understood from Figure 2, which plots the survival probability 
of electron neutrinos for normal and inverted mass hierarchies. 

5.2.1 Small and Icirge mixing angle approximations 

We can first distinguish four limiting cases depending on the type of mass hierarchy 
and the value of the 6'i3 mixing angle (sec equations 5 and 6 and Figure 2): 

1. n.h.-L for normal mass hierarchy and large ^13, where Pee ~ and Pee ~ 0.7, 

2. n.h.-S for normal hierarchy and small ^13, where Pee ~ 0.3 and Pee ~ 0.7, 

3. i.h.-L for inverted mass hierarchy and large ^13, where Pee ~ 0.3 and Pee ~ 0, 

4. i.h.-S for inverted mass hierarchy and small ^13, where Pee ~ 0.3 and Pee ~ 0.7, 

and where "small" (resp. "large" ) mean sin ^13 < 2 X 10-6 (resp. sin^ ^13 > 3 x 10"^). 
In the intermediate mixing angle region, there survival probability depends on the 
neutrino energy. Note that the no oscillation case corresponds to Pee = 1 and Pee = 1. 

We can therefore summarize the situation as follows: the suppression of the i/e flavor 
is maximal in the case of a normal mass hierarchy and large ^13 mixing angle. This 
corresponds to the total conversion of 1/^8 into ^';u,r's as they traverse the supernova 



14 



x10 



s1200 




0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 

time (s) 




time < 0.47 s 



10 15 20 25 30 35 40 45 50 
Energy (MeV) 




5 10 15 20 25 30 35 40 45 50 

Energy (MeV) 



Figure 5: Illustration of the prediction of [16] without neutrino flavor oscillations: 
prediction of the neutrino flux at Earth at a distance of 10 kpc (top), time-integrated 
energy spectrum corresponding to the first 270 ms after the bounce (bottom-left) and 
luminosity spectrum of Vf. (solid line), (dashed line) and v^^r-, ^ii,T (dotted line) 
neutrinos at 100 ms after the shock breakout (bottom-right). We treat the times 
before t=240 ms as "burst" and after as "cooling". 
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Burst (t < 240 ms) - 3 kton detector 

Reaction No Oscillation (n.h.) Oscillation (i.h.) 

oscillation Large ^13 Small ^13 Large ^13 Small ^13 

Elastic 







2 


<1 


1 


1 


1 




Pee" 


<1 


<1 


<1 


<1 


<1 






<1 


<1 


<1 


<1 


<1 




(P^ + Pr) e" 


<1 


<1 


<1 


<1 


<1 




total V e~ 


2 


<1 


1 


1 


1 


Absorption 














CC 


i/e ^°Ar 


8 


2 


4 


4 


4 




Pe 4°Ar 


<1 


<1 


<1 


<1 


<1 


NC 


z/4"Ar 


1 


1 


1 


1 


1 




p4°Ar 


<1 


<1 


<1 


<1 


<1 


Total 




11 


3 


6 


6 


6 



Tabic 2: Expected neutrino events in a 3 kton detector during the burst (see text for 
definition). Neutrino oscillations with matter effects are included for different mass 
hierarchies and 6*13 values. 



matter. For the other cases, the oscillation reduces the expected flux by a factor 3 
compared to the non-oscillation case. 

Tables 2 and 3 show the expected number of events during the burst in a 3 kton 
and 70 kton hquid Argon TPC detector, respectively, for the four oscillation scenarios 
and compares them to the non oscillation case. For completeness we have separated 
the contributions from z/g's, f^^r^ and similarly for antineutrinos. However, it should 
be recalled that four channels are in fact accessible experimentally: the Uf. CC, the Pe 
CC, the NC and the elastic scattering. 

With the assumed distance of 10 kpc, the burst produces in a 3 kton detector and 
in absence of neutrino oscillations 2 elastic events, 8 CC and 1 NC event for a total 
of 11 events. In case of the 70 kton detector, the rates are 37 elastic events, 192 CC 
and 39 NC event for a total of 268 events. 

We consider first the elastic events. As already mentioned, oscillation has a strong 
effect on the z/g burst. The z/g will be suppressed by a factor 3 for three of the four 
oscillation scenarios we have considered, or can even be maximally suppressed in the 
case of normal mass hierarchy and large ^13. Neutrino flavor oscillation is of course 
a unitary process, hence, the oscillated will be transformed into other flavors. 
However, since the cross-section for the elastic z/g process is higher than for other 
processes, we expect a reduction of the number of elastic events. Indeed, as visible in 
Tables 2 and 3, the elastic burst signal is highly suppressed for the oscillation scenarios. 
With our deflnition of burst, with the assumption of a 40 ms duration, there is in fact 
a smaller component of i/^u's, v^s, P^^'s, P,-'s and Pg's (see Figure 5 ). Finally, the elastic 
events are suppressed by a factor roughly 4 for the case of a normal mass hierarchy 
and large ^13 mixing angle, and roughly by a factor 2 for the other cases. 
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Burst (t < 240 ms) - 70 kton detector 
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86 




Pe ^°Ar 


<1 


<1 
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1 


<1 


NC 


z/4<JAr 


28 


28 


28 


28 


28 




p40Ar 


n 


11 


11 


11 


11 


Total 




268 


91 


144 


145 


144 



Tabic 3: Expected neutrino events in a 70 kton detector during the burst (see text 
for definition). Neutrino oscillations with matter effects are included for different mass 
hierarchies and 6*13 values. 



We now turn to the nuclear charged current events. The z^e CC reaction is in 
principle a golden channel to study the burst. Because of the oscillation of into 
other fiavors (not completely compensated by oscillation of other flavors into v^. the 

size of the peak is reduced compared to the non-oscillation case: the suppression is 
almost a factor 5 for the case of a normal mass hierarchy and large ^13 mixing angle, 
and roughly a factor 2.5 for the other cases. 

Finally, we can consider the neutral current events. Obviously they are not affected 
by oscillations, since neutral currents are flavor symmetric. They hence constitute an 
excellent probe for the supernova properties independent of the neutrino oscillation 
intrinsic properties. We therefore stress the importance of this signal to decouple 
supernova from oscillation physics. Unfortunately the rates in the 3 kton detector 
are rather modest: at the level of one event (independent of the oscillation scenario). 
For the case of a 70 kton detector, we expect in total 39 events (independent of the 
oscillation scenario). 

The information about the shock breakout mechanism that can be extracted from 
the 3 kton detector is hence marginal, in particular in the case of the reduced number 
of expected events because of oscillation effects. A supernova could explode at a closer 
distance, however, this is an unlikely event. In contrast, a 70 kton detector could 
provide hundreds of events, enough statistics to study the neutrinos from the burst. 

5.2.2 Study of energy spectra 

The importance of the statistics can be appreciated when one studies the energy spectra 
of the events. Sufficient statistics would allow to study the energy spectra of the events 
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which provides an extra signature to separate supernova from oscillation physics. 

Figures 6-9 contain the time evolution of the expected events in the burst for 
the different detection channels, namely i^e CC, Pg CC, elastic and neutral currents, 
normalized to 70 kton detector and the corresponding time integrated event spectra. 
The non oscillation and the four different oscillation cases are compared in the plots. 
For the z/g CC, CC and clastic events we also show the event energy distributions. 
For the neutral current events, the energy cannot be determined due to the neutrino 
in the final state. 

The reduction of the peak due to the effect of oscillations can be seen in figure 
6. The suppression is maximal for normal hierarchy and large ^13 mixing angle due 
to the total conversion of into f^^,-. At the same time, the energy spectrum moves 
slightly to higher neutrino energy values. 

Elastic events (figure 8) are also sensitive to the oscillations with a reduction of the 
shock breakout peak. However, the energy spectrum remains almost unchanged. 

The neutral current events (figure 9) are not affected by neutrino oscillations and 
the time evolution clearly shows the shock breakout peak coming from z/g events. 

Overall, we note that the energy distribution of the burst is marginally affected by 
oscillations. 

5.2.3 Consequence of the neutrino burst detection 

We can hence conclude that the burst signal is only affected in rate by the effect of 
neutrino oscillation and the effect on the shape is marginal. This means that: 

• a clear SN burst peak should be visible in the neutral current sample, regardless of 
the neutrino oscillation physics. This signature should provide a clear indication 
for the nature of the supernova explosion. The size of the peak determines the 
burst absolute neutrino flux from the supernova. 

• the oscillation physics occurring during the early phase of the supernova can be 
studied by detection of the and Pg CC channels. 

• the number of observed v^, CC observed in the burst can distinguish between 
the non oscillation case, the n.h.-L oscillation and the three other oscillation 
scenarios. The n.h.-S, i.h.-L and i.h.-S cases cannot be distinguished by Ue CC. 

• a roughly factor 4 enhancement of observed Pg CC observed in the burst would 
be a clear indication for an oscillation matter enhancement only possible in the 
inverted mass hierarchy with large mixing angle. Otherwise, a roughly factor 2 
enhancement is expected from the effect of oscillations (see Figure 7) . It should 
be stressed that this effect is not visible even for a 70 kton detector. 

• elastic events can help distinguish between the non oscillation case, the n.h.-L 
oscillation case and the three other oscillation scenario. 
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Figure 6: SN burst: time evolution of the z/g CC event rate (top) and the corresponding 
time integrated event spectra (bottom). The non oscillation (dark curve), n.h.-L (light 
curve) and n.h.-S, i.h.-L and i.h.-S (dark region) cases are considered. 
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Figure 7: SN burst: time evolution of the z/g CC event rate (top) and the corresponding 
time integrated event spectra (bottom). The non oscillation (dark curve), i.h.-L (light 
curve) and n.h.-L, n.h.-S and i.h.-S (gray region) cases are considered. 
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Figure 8: SN burst: time evolution of the elastic event rate (top) and the corresponding 
time integrated event spectra (bottom). The non oscillation (dark curve), n.h.-L (light 
curve) and n.h.-S, i.h.-L and i.h.-S (dark region) cases are considered. 
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Figure 9: SN burst: time evolution of the NC event rate. 

5.2.4 General mixing case — intermediate mixing angle 

For the intermediate mixing angle region defined approximately by ~ 10~^ < sin^ ^13 < 
2 X 10~^, the number of expected events varies strongly with the mixing angle. In fact, 
in this region, we expect that the neutrino oscillation effect depends on the neutrino 
energy within the range up to 100 McV relevant for supernova neutrinos. 

The actual variation of the number of events as a function of the ^^13 value arc 
summarized in Figure 10 for the four detection channels. Solid lines correspond to 
the normal mass hierarchy case, dotted lines to the inverted mass hierarchy case and 
dashed lines to the non oscillation case. We consider events in burst and the plots are 
normalized to a 70 kton detector. 

This variation can be exploited to measure the mixing angle if it lies in this inter- 
mediate region. Quantitative estimate of the precision with which the angle can be 
determined will be described elsewhere [23]. 

5.3 Comparison with other experiments 

All experiments will be affected in the same way by the effect of neutrino fiavor oscil- 
lation in their expected rates of Ve events. In order to compare experiments, we hence 
restrict ourself to the case where neutrino oscillations are neglected and compare rates 
in Superkamiokande and SNO with our results on Argon. Results can be readily scaled 
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Figure 10: Expected number of events in a 70 kton detector during the burst (see text 
for definition) as a function of sin^ ^13. The different neutrino interaction processes are 
plotted separately. Solid lines correspond to the n.h. oscillation case, dotted lines to 
the i.h. case and dashed lines to the non oscillation case. 
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to include neutrino oscillation effects. 

The integrated number of z/g events as a function of time for elastic and absorption 
processes without including oscillation effects are plotted in figure 11. The main con- 
tribution comes from the CC channel. Figure 12 shows for comparison similar plots 
for SuperKamiokande (SK) and SNO detectors. 

The SK volume for supernova neutrino detection is 32 kton of light water. The 
dominant reaction is the charged-current absorption of z/g neutrinos on free protons 
{ueP ne+). Neutrino scattering on electrons (i/e~ — > i/e~) also contributes for all 
neutrino species. For the study of the shock breakout we calculate the event rates in 
SK due to only the elastic scattering. For this process we use the neutrino electron 
scattering cross section given in [24] with a threshold of ~ 5 MeV. If there is large 
mixing angle between 1/^ and or i/^ neutrinos, the early i/g signal will be reduced 
significantly in SK. 

The SNO detector is a light and heavy water neutrino detector containing 1 kton 
of D2O surrounded by a cavity of 1.6 kton of light water. The main reactions for 
i/e neutrinos are the elastic scattering and the NC and CC reactions on deuterium 
{ued npue and Ued — > ppe"). The cross sections used to compute these events are 
taken from [25]. 

Figure 13 shows the comparison between the expected number of events from the z/g 
burst in the three experiments. A total of 10 events arc expected in the 3kton ICARUS 
in the 40 ms after the bounce: 8 events from the absorption channel and 2 from the 
elastic channel. The contribution of SK and SNO for the same interval of time is 20 
and 7 events, respectively. 

Summarizing, we observe that the event statistics that would be obtained in Su- 
perKamiokande, SNO and the 3 kton ICARUS arc similar. In this case, SNO and 
ICARUS benefit largely from their direct sensitivity to z/g via charged current interac- 
tions on resp. deuteron and Argon, while in SuperKamiokande one has to rely on elastic 
scattering on electrons, which has a much smaller cross-section. Hence, smaller detec- 
tors hke SNO and ICARUS can compare favorably with the largest SuperKamiokande 
in the study of the breakout phase. 

We also conclude that given that neutrino oscillations is an established fact [2], it 
will be difficult to study the neutrino burst within the current round of experiments 
given the reduction of statistics following the conversion of z/g into other fiavors. In the 
context of the liquid Argon TPC, it appears that a 70 kton detector will be mandatory 
in order to collect sufficient events to study the neutrino burst. 

6 Cooling phase 

In this section we will refer to the main neutrino emission phase in the core collapse 
supernovae: the cool-off of the star. 

As seen in previous sections, the energy spectra of the expected neutrinos at the 
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Figure 11: Integrated number of z/g events as a function of time for the elastic and 
absorption channels. The thick solid line corresponds to the total number of events for 
a 3 kton detector. No oscillation effects are included. 
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Figure 12: Integrated number of events as a function of time expected in the SK 
(left) and SNO (right) detectors. No oscillation effects are included. 
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Figure 13: Comparison of the expected number of events from the z/g breakout burst 
in a 3 kton ICARUS, SK and SNO experiments. In the first 40 ms after bounce a total 
of 20 events are expected in SK, 10 in ICARUS and 7 in SNO. No oscillation effects 
are included. 
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detector is modified by oscillations among flavors since the average neutrino energies 
are flavor dependent and the interaction cross sections are different for different flavors. 
The variations of the neutrino energy will depend on the adiabaticity conditions on the 
H and L resonances in the supernova and on the neutrino oscillation parameters. Four 
different cases are considered in this analysis depending on the type of mass hierarchy 
and the value of the 9 is mixing angle. 

6.1 Hierarchical energy distribution — scenario I 

We first consider the prediction from Ref [17] which we parameterize as Fermi-Dirac 
distributions with {E^^) = 11 MeV, {EpJ = 16 MeV, (^^^^) = (Ep^^) = 25 MeV and 
luminosity equipartition among flavors. In this kind of scenario where the energies are 
rather different (hierarchical) between flavors, the effect of the oscillation is the most 
pronounced. Indeed, oscillations will as expected harden the spectrum, which has a 
big incidence of the Ve CC channel. 

Table 4 contains the number of neutrino events from a supernova at 10 kpc ex- 
pected in the 3 kton detector. Elastic, CC and NC processes have been considered 
independently and oscillation and non oscillation cases have been computed for normal 
and inverted hierarchies. However, we recall that four channels are in fact accessible 
experimentally: the Ue CC, the Ue CC, the NC and the elastic scattering. 

We note that: 

• the number of observed CC observed can distinguish between (a) the non 
oscillation case, (b) the n.h.-L case and (c) the three other oscillation scenarios. 
n.h.-S, i.h.-L and i.h.-S cannot be distinguished by CC. 

• a roughly factor 4 enhancement of observed Pg CC observed would be a clear in- 
dication for an oscillation matter enhancement only possible in the inverted mass 
hierarchy with large mixing angle. Otherwise, a roughly factor 2 enhancement is 
expected from the effect of oscillations. 

As was mentioned in section 3, normal and inverted hierarchies are degenerated for 
small in all the cases. In the neutrino (antineutrino) channel and i.h. (n.h.) the 
conversion is adiabatic and results are independent on ^13. Maximal conversion occurs 
for large mixing angle in the neutrino (antineutrino) channel for n.h. (i.h.). 

We stress again that the flux and flux have different information about the 
neutrino oscillation parameters. The i.h.-L and i.h.-S are distinguishable from Pg events 
but are not from i/g events. Then, a good separation of these events is important in 
order to distinguish between models. 

Figures 14, 15 and 16 show these effects in the neutrino energy spectra. The plots 
on the top correspond to the expected number of neutrinos arriving at Earth and the 
bottom plots are the corresponding number of events detected in a 3 kton detector. 
The four oscillation cases are considered in the distributions. 
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Scenario I 



Reaction 




No 


Oscillation (n.h.) 


Oscillation (i.h.) 






oscillation 


Large 6113 
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Large 6^13 


Small 6^13 


Elastic 


















20 


20 


20 


20 


20 




z/g e~ 
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8 




(z/^ + Ur) e" 
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7 


7 




{Uf, + v^) e" 
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6 


6 




total z/e~ 


41 


41 


41 


41 


41 


Absorption 














CC 


z/e 40 Ar 


188 


962 


730 


730 


730 




z/e 40 Ar 


15 


33 


33 


75 


33 


NC 


z/ 40 Ar 


492 


492 


492 


492 


492 




z/40Ar 


419 


419 


419 


419 


419 


Total 




1155 


1947 


1715 


1757 


1715 



Table 4: Expected neutrino events in a 3 kton detector including neutrino oscillations 
with matter effects inside the supernova. The flux parameters are {E^^) = 11 MeV, 
(-EpJ = 16 MeV, (-Ey^^) = (-Ep^^) = 25 MeV and luminosity equipartition among 
flavors is assumed. 
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Figure 14: Scenario I: Number of z/g neutrinos arriving at Earth (top) and expected 
number of events in the 3 kton detector (bottom) for z/g neutrinos. The non oscillation 
and the four oscillation cases have been taken into account. The CC interaction process 
(left) and elastic scattering (right) are split in the figure. 
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Figure 15: Scenario I: Number of z/g neutrinos arriving at Earth (top) and expected 
number of events in the 3 kton detector (bottom) for neutrinos. The non oscillation 
and the four oscillation cases have been taken into account. The CC interaction process 
(left) and elastic scattering (right) are split in the figure. 
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Figure 16: Scenario I: Number of neutrinos arriving at Earth (top) and expected num- 
ber of events in the 3 kton detector (bottom) for i/^^,- (left) and v^^r (right) neutrinos. 
The non oscillation and the four oscillation cases have been taken into account. 
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Due to the total conversion z/^^^- — > v^. for the n.h.-L case, the v^. energy spectra 
is harder and this leads a huge increase of the expected events due to the quadratic 
dependence of the CC cross section with energy. The same effect can be seen for v^. 
events and i.h.-L case. The elastic processes are less sensitive to the oscillations and 
smaller number of events is expected. Nevertheless, the energy spectrum is modified 
specially for v^^^ events, moving to lower values of energy due to the neutrino mixing. 

6.1.1 General results as a function of the mixing angle 

The explicit variation of the expected rates and energy spectra with the ^13 angle can be 
computed considering the dependence of the jump probability with this angle (section 
3). 

Figure 17 shows the the variation of the neutrino rates as a function of sin^ 6'i3. We 
compare the expected results for the cases of n.h. (solid line), i.h. (dotted line) and 
no oscillation (dashed hue) . CC events are very sensitive to the change on the energy 
spectra due to oscillations. A clear increase on the number of events is expected, being 
of a factor (4-5) for CC events and a factor 2 for z/g CC interactions. The main 
variations with the angle ^13 are expected for the CC channel in the case of normal 
hierarchy and for the Pg CC channel for inverted hierarchy. The number of CC 
events increases 30% from small to large ^13 values and Pg CC events are doubled. 
Small variations are expected for the rest of channels. 

Although the number of events does not change for some reactions, the energy 
spectra are affected by oscillations and in particular by the value of the ^13 angle. 

6.2 Other cooling scenarios — scenarios II— V 

The signature of neutrino oscillations depends on the flavor-dependent primary super- 
nova neutrino fluxes. Recent SN simulations [7, 18] found that the differences between 
the average energies of Vfx.Ti J^ij.,t and Pe are in the range 0-20%, with 10% being a typ- 
ical number. The himinosities of all neutrino flavors are expected to be approximately 
equal within a factor two. 

In order to consider the range of possibilities for the source fluxes we compute the 
expected events for different values of the neutrino average energies and luminosities, as 
summarized in table 1. The results from the different scenarios are contained in tables 
5-8. Table 5 is obtained using the flux parameters found by the Livermore group [10] . 

In table 6 we take the results from [7, 18] and equipartition luminosities. Finally, 
tables 7 and 8 show the expected rates considering the same average energies than used 
in table 6 but for two extreme cases of luminosity partition: L^^ — Lp^ = 2 Lj^^ (table 
7) and L^^ = Lp^ = 0.5 L^^ (table 8), with i/^ each of z/^, Ur, P/i and neutrino flavors. 

Clearly, scenarios III-V distinguish themselves from the scenarios I-II in their pre- 
diction of a quite flavor degenerate spectra (see figures 14, 15 and 16). In order to 
illustrate this, we show in Figure 18 how these affect the neutrino energy spectra. The 
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Figure 17: Expected number of events in a 3 kton detector as a function of sin^^is. 
The different neutrino interaction processes expected are plotted separately. Solid lines 
correspond to the n.h. oscillation case, dotted lines to the i.h. case and dashed lines 
to the non oscillation case. 
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total u e~ 


40 


40 


40 


40 


40 


Absorption 














CC 


Vg 40 Ar 


264 


805 


644 


644 


644 




Pe ^°Ar 


16 


28 


28 


57 


28 


NC 


z/4<JAr 


415 


415 


415 


415 


415 




p40Ar 


355 


355 


355 


355 


355 


Total 




1090 


1643 


1482 


1511 


1482 



Tabic 5: Expected neutrino events in a 3 kton detector. The flux parameters are i^E^,^ 
= 13 McV, {E^J = 16 MeV, {E^^J = {Ep^J = 23 MeV and luminosity equipartition 
among flavors is assumed. 



upper plot corresponds to the expected number of neutrinos arriving at Earth and the 
bottom plot is the event energy spectrum normahzed to a 3 kton detector. The three 
scenarios III-V are considered in the distributions, each time for the non oscillation 
and for the oscillation case with a normal mass hierarchy and large mixing angle. 

6.3 Discrimination between the various scenarios 

The comparison between the different scenarios is summarized in figure 19. The ex- 
pected number of events in the four neutrino detection channels are plotted for the five 
scenarios of the SN neutrino fiux parameters considered in this study. Moreover, the 
non oscillation and the four oscillation cases are plotted for every scenario. The plots 
are normalized to a supernova at 10 kpc and a 3 kton detector. The errors indicated 
in the picture correspond to statistical errors and give a visual feeling on the ability to 
separate the various supernova (I~V) and oscillation scenarios (no osc, n.h.-L, n.h.-S, 
i.h.-L, i.h.-S). 

NC events are sensitive to the primary source parameters and are independent 
to the oscillation effects. They give a direct account for the total luminosity of the 
supernova and of the hardness of the neutrino spectra. Hence, they provide the means 
to discriminate between the various supernova physics scenarios if the binding energy 
is assumed to be known. This is illustrated in the bottom plot of Figure 19. The rates 
of the nuclear processes are such that excellent statistical precision is achieved. They 
amount to 911 ±30, 770 ±27, 422 ±20, 387 ±20, 450 ±21 events for resp. the scenarios 
I-V, where the error quoted is the statistical error on the number of events. 
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Scenario III 



Reaction 




No 


Oscillation (n.h.) 


Oscillation (i.h.) 






oscillation 


Large ^13 


Small 


Large ^13 


Small ^13 


Elcistic 
















i/g e 


20 


20 


20 


20 


20 






8 


8 


8 


8 


8 




(i/^ + I'r) e" 


7 


7 


7 


7 


7 






5 


5 


5 


5 


5 




total u e" 


40 


40 


40 


40 


40 


Absorption 














CC 




264 


482 


417 


417 


417 




z/e ^°Ar 


16 


18 


18 


23 


18 


NC 


z/*"Ar 


215 


215 


215 


215 


215 




i"Ar 


207 


207 


207 


207 


207 


Total 




742 


962 


897 


902 


897 



Table 6: Expected neutrino events in a 3 kton detector. The flux parameters are {Ey^) 
= 13 MeV, (£;pj = 16 MeV, {E^^^,) = (-Ep^,,) = 17.6 MeV and luminosity equipartition 
among flavors is assumed. 









Scenario IV 






Reaction 




No 


Oscillat 


ion (n.h.) 


Oscillat 


ion (i.h.) 






oscillation 


Large 6113 


Small 6*13 


Large 6113 


Small 6*13 


Elastic 


















29 


15 


19 


19 


19 






12 


10 


10 


6 


10 




(z/^ + Vr) e" 


5 


8 


7 


7 


7 




(z/^ + Vr) e" 


4 


5 


5 


6 


5 




total v e~ 


50 


38 


41 


38 


41 


Absorption 














CC 


z/e ^°Ar 


397 


362 


372 


372 


372 




z/e ^°Ar 


24 


22 


22 


17 


22 


NC 


z/4"Ar 


188 


188 


188 


188 


188 






199 


199 


199 


199 


199 


Total 




858 


809 


822 


814 


822 



Table 7: Expected neutrino events in a 3 kton detector. The flux parameters are {Ey^) 
= 13 MeV, {Ej,^) = 16 MeV, {E^^J = {Ep^J = 17.6 MeV and luminosities = Lp^ 
— 2 with i/x each of 1/^, 1/^, and D-r neutrino flavors are considered. 
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No osc (seen III) 
n.h.-L (seen III) 
No ose (seen IV) 
n.h.-L (seen IV) 
No ose (seen V) 
n.h.-L (seen V) 



60 80 100 

energy (MeV) 




No ose (seen III) 
n.h.-L (seen III) 
IV) 
IV) 
V) 
V) 



100 



energy (MeV) 



Figure 18: Scenarios III-V: number of neutrinos arriving at Earth (top) and expected 
number of events in the 3 kton detector (bottom) for z/g neutrinos. The three scenarios 
III-V are considered in the distributions, each time for the non oscillation and n.h.-L 
oscillation cases. 
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Scenario V 



l-e /^T1 /'"ATI 




iNO 




1011 ^^Il.Il.J 




lOll ^^l.Il.J 






OfeCllldtlOIl 


Liarge ^13 


Ollldll t/13 


Liarge fia 


Ollldll t/13 


Jiilastic 
















Ug e" 


1 
IZ 


ZO 


on 
zu 


on 
zu 


on 
zu 






r 



n 
D 


D 


1 n 


D 






Q 

o 


L) 


7 


7 


7 




/ — 1 — \ — 

{Uf, + u^) e 


7 


6 


6 


5 


6 




total u e~ 


32 


41 


39 


42 


39 


Absorption 














CC 


Vg 40 Ar 


159 


578 


453 


453 


453 




Pe ^°Ar 


10 


15 


15 


28 


15 


NC 


z/4<JAr 


236 


236 


236 


236 


236 




p40Ar 


214 


214 


214 


214 


214 


Total 




651 


1084 


957 


973 


957 



Tabic 8: Expected neutrino events in a 3 kton detector. The flux parameters are {E^^ 
= 13 MeV, (EpJ = 16 MeV, (^^^,,) = (^p^.,) = 17.6 MeV and luminosities L^,^ = Lp^ 
= 0.5 \jy^ with t/j. each of t/^, i^t-, P^ and Vt neutrino flavors are considered. 



The supernova and neutrino oscillation physics can then be studied by the simul- 
taneous observation of NC, v^. CC, Pg CC and elastic channels. We observe: 

• In the hierarchical scenarios (1-11) the enhancement of due to oscillation is 
highly signiflcant. As already mentioned, this is the result of the hardening of 
the spectrum in all cases of oscillations. 

— The effect is of course the strongest in the case of n.h.-L for the ^-nd in 
the case of i.h.-L for Pg. 

— The normal and inverted hierarchy yield the same number of events in the 
case of small mixing. 

• In the non-hierarchical scenarios (III-V) the enhancement of Vf. due to oscillation 
is signiflcant compared to the non-oscillation case, except for scenario (IV) where 
the quasi-degenerate energies and the choice of luminosities Lj^^ = Lp^ = 2 
with each of z/^, v^^ P^ and neutrino flavors just washes out the oscillation 
(the disappearing v^, — > Vj. are compensated by the appearing Vj. — > v^. 

• in the scenario IV this ambiguity in the number of events can be solved by a 
study of the event energy distribution (see Figure 18). Indeed, the flgure clearly 
shows that the spectrum is indeed harder for the oscillation case but that the 
total number of events is roughly the same. 
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Figure 19: Expected number of events: the four neutrino processes are considered for 
different SN scenarios and oscillation conditions (see text). The statistical errors are 
included in the picture. 
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7 Conclusions 



An analysis of the neutrino events expected in a liquid Argon TPC (ICARUS-like) 
from a future galactic supernova has been performed showing the high capabilities of 
the detector to observe neutrinos from stellar collapses and extract information about 
their properties. 

Four detection channels have been considered: elastic scattering on atomic electrons 
from all neutrino species, v^, charged current absorption on Ar with production of 
excited Pg charged current absorption on Ar with production of excited CI and 
neutral current interactions on Ar from all neutrino flavors. 

We give a detailed account of the dependence with the mixing angle ^13 and the type 
of mass hierarchy. Neutrino oscillations change significantly the expected supernova 
neutrino rates and energy spectra. The matter effects inside the supernova must be 
included to predict how the spectra of the neutrino burst is modified by oscillations. 

Information about the supernova explosion mechanism can be obtained from the 
study of the neutrino burst. The early phase neutrino signal on Argon has been com- 
puted using the fluxes from a realistic core collapse supernova model [16], and including 
oscillation effects. Thanks to its clean identification and high sensitivity to neutrinos, 
a liquid Argon TPC can provide fundamental information about the shock breakout 
mechanism. We find that a 70 kton detector is best suited to study the SN burst, given 
its statistics even in case of the z/g suppression expected in case of oscillations. 

For the cooling phase of the supernova, we have adopted a simple Fermi-Dirac 
parameterization of the neutrino fluxes. We selected six different scenarios in order 
to take into account the uncertainties on the knowledge of the supernova physics and 
estimate their effect on our analysis. For the study of the cooling phase, a 3 kton 
detector is sufficient to collect hundreds of events which could provide fundamental 
information on supernova and neutrino oscillation physics. 

We discussed how the various scenarios can be distinguished. In all cases, we find 
that the simultaneous observation of the four channels, in particular nuclear charged 
and neutral current processes, is fundamental in order to decouple the supernova from 
the neutrino oscillation physics. 

The observed event rates and the event energy spectra of the four classes of events 
form a set of observables sensitive to the supernova parameters and to the intrinsic neu- 
trino oscillation parameters (^13 and mass hierarchy). We have qualitatively illustrated 
this with several examples. 

In order to understand quantitatively the sensitivity to such parameters, one should 
study the classification of actual detected events into the four channels and perform a 
global fit to the observed distributions. Such an analysis will be the subject of a future 
work [23]. 
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